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O i Abstract: We have developed a physical model for the non-perturbative x- 

shape of parton density functions in the proton, based on Gaussian fluctuations 

in momenta, and quantum fluctuations of the proton into meson-baryon pairs. 

^ . The model describes the proton structure function and gives a natural explana- 

^ ! tion of observed quark asymmetries, such as the difference between the anti-up 

^ I and anti-down sea quark distributions and between the up and down valence 

^ I distributions. We find an asymmetry in the momentum distribution of strange 

^ ■ and anti-strange quarks in the nucleon, large enough to reduce the NuTeV 

^ ■ anomaly to a level which does not give a significant indication of physics be- 

^ . yond the standard model. We also consider charmed fluctuations, and show 

OO I that they can explain the excess at large x in the EMC F| data. 
O . 
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The low-scale parton density functions give a description of the hadron at a non- 
^ ■ perturbative level. The conventional approach to these functions is to make parameteriza- 
>■ ! tions using some more or less arbitrary functional forms, based on data from deep inelastic 
' scattering and hadron collision experiments. Another approach, however, is to start from 
^ I some ideas of the behavior of partons in the non-perturbative hadron, and build a model 
based on that behavior. The advantage with this approach is that the successes and fail- 
ures of such a model allows us to get insight into the non-perturbative QCD dynamics. 
The model presented here, and described in detail in [H |21, describes the F2 structure 
function of the proton, as well as sea quark asymmetries of the nucleon. Most noteworthy, 
our model predicts an asymmetry between the momentum distributions of strange and 
anti-strange quarks in the nucleon of the same order as the newly reported results from 
NuTeV 0. The model also suggests an intrinsic charm component in the proton. 

This work extends the model previously presented in The model gives the four- 
momentum k of a single probed valence parton (see Fig. ^ for definitions of momenta) by 
assuming that, in the nucleon rest frame, the shape of the momentum distribution for a 
parton of type i and mass mj can be taken as a Gaussian 

f,{k) = N{a,,m,) exp {- [{k^ - nnf + kl + k'y + k^] /2a^} (1) 



^This is an extended version of the talk "Quark Asymmetries in Nuclcons", given at the XIII Interna- 
tional Workshop on Deep Inelastic Scattering, Madison, USA, April 27-May 1, 2005 
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FIGURE 1. Illustration of the processes probing (a) a valence parton in the proton and (b) a 
sea parton in a hadronic fluctuation (letters are four- momenta) . (c) shows the resulting parton 
distributions at the starting scale Qq. 

which may be motivated as a result of the many interactions binding the parton in the nu- 
cleon. The width of the distribution should be of order hundred MeV from the Heisenberg 
uncertainty relation applied to the nucleon size, i.e. cTi = l/d^- The momentum fraction x 
of the parton is then defined as the light-cone fraction x = k^/p^. We impose constraints 
on the final-state momenta in order to obtain a kinematically allowed final state, which 
also ensures that < x < 1 and fi{x) — for x — 1. Using a Monte Carlo method these 
parton distributions are integrated numerically without approximations. 

To describe the dynamics of the sea partons, we note that the appropriate basis for the 
non-perturbative dynamics of the bound state nucleon should be hadronic. Therefore we 
consider hadronic fluctuations, for the proton 



Probing a parton i in a hadron if of a baryon- meson fluctuation \BM) (see Fig-QJ^) gives a 
sea parton with light-cone fraction x = xh Xi of the target proton. The momentum of the 
probed hadron is given by a similar Gaussian, but with a separate width parameter cth- 
Also here, kinematic constraints ensure that we get a physically allowed final state. The 
procedure gives xh ~ Mh/ (Mb+Mm), i-e. the heavier baryon gets a harder spectrum than 
the lighter meson. The normalization of the sea distributions is given by the normalization 
coefficients a%M ^'i- ©• These cannot be calculated from first principles in QCD and 
are therefore taken as free parameters to be fitted using experimental data. 

The resulting valence and sea parton x-distributions apply at a low scale Qq, and 
the distributions at higher are obtained using perturbative QCD evolution at next-to- 
leading order. 

The model has in total four shape parameters and three normalization parameters, plus 
the starting scale, to determine the parton densities u, d, g, u, d, s, s. These are (with 



\p) = aobo) + ttpTrobTr") + a„^+|n7r+) + . . . + aAK\AK+) + ... 



(2) 
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FIGURE 2. The proton structure function F2{x,Q'^) for large x values; NMC and BCDMS 
data El compared to our model, also showing the results of ±20% variations of the width 
parameters cr^ and ad for the u and d valence distributions. 
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The resulting parton densities are shown in Fig. ^c). 

In order to fix the values of the model parameters, we make a global fit using several 
experimental data sets: Fixed-target F2 data to fix large-x (valence) distributions (Fig.|21); 
HERA F2 data for the gluon distribution width and the starting scale Qo! cF/u- asymmetry- 
data for the normalizations of the \p'K^) and |?T.7r"'") fluctuations (see Fig. 4); and strange 
sea data to fix the normalization of fluctuations including strange quarks (see Fig. 5a). 
We have also compared with charge asymmetry data as a cross-check on the ratio of 
Gaussian widths for the u and d valence quark distributions (Fig. 3). It is interesting to 
note that this simple model can describe such a wealth of different data with just one or 
two parameters per data set. 




FIGURE 3. The charge asymmetry 
for leptons from VF^-decays in pp col- 
lisions at the Tevatron ijj compared 
to our model, with best-fit parameters 
and a 20% reduced width of the valence 
d quark distribution. 
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FIGURE 4. Comparison 
between our model and data 
from the E866/NuSea collab- 
oration 0: (a) u{x)/d{x) (b) 
xd{x) — xu{x). The full line 
uses the physical pion mass, 
while the dashed line uses an 
effective pions mass rrf^ = 
400 MeV as discussed in the 
text. 



In our model, the shape difference between the valence u and d distributions in the 
proton, apparent from the charge asymmetry data, is described as different Gaussian 
widths. This would correspond to a larger effective volume in the proton for d quarks than 
for u quarks, an effect which could conceivably be explained by Pauli blocking of the u 
quarks. 

Since the proton can fluctuate to 7r° and 7r+ by |p7r°) and InTT"*"), but to vr^ only by 
the heavier IA+^tt"), we get an excess of d over u in the proton sea. Interestingly, the fit 
to data improves when we use a larger effective pion mass of 400 MeV (see Fig. 4). This 
might indicate that we have a surprisingly large coupling to heavier p mesons, or that one 
should use a more generic meson mass rather than the very light pion. 

The lightest strange fluctuation is \KK^). If we let this implicitly include also heavier 
strange meson-baryon fluctuations, we can fit the normalization to strange sea data 
(see Fig. 5a). The result corresponds to ^^{xs + xs)dx/ {xu + xd)dx ~ 0.5, in agreement 
with standard parton density parameterizations. We note that this indicates a normal- 
ization oc 1/AMbm = I /{Mb + Mm — Mp) rather than oc as expected from 
old-fashioned perturbation theory. The fluctuation parameters are taken from the light sea 
results, an = 100 MeV and aq = a^"^"^ = 170 MeV as discussed in p. 

Since the s quark is in the heavier baryon A and the s quark is in the lighter meson , 
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FIGURE 5. (a) CCFR data P on the strange sea distribution {xs{x) + xs{x))/2 compared 
to our model based on \AK) fluctuations with different normalizations, (b) The strange sea 
asymmetry s~{x) = xs{x) — xs{x) (at = 20 GeV^) from the model and combined with the 
function F{x) accounting for NuTeV's analysis giving A sin = dxs^{x)F{x) = -0.0017. 
The uncertainty bands correspond to the uncertainties for S~ and Asin^^vi^ quoted in the text. 
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FIGURE 6. (a) Comparison between the 
strange and charm sea obtained from our model 
with the inclusion of the h.cD fluctuation. The 
normalization is here taken to be oc 1/{Mb + 
Mm — Mp) as suggested by strange sea data, (b) 
The charm quark distributions from the p J/ip 
fluctuation (the c and c distributions are identi- 
cal with this fluctuation). Note that the distri- 
butions are very similar to those from the AqD 
7 distribution, except for the small asymmetry be- 
X tween c and c. 

we get a non-zero asymmetry = dx[xs{x) — xs{x)] in the momentum distribution of 
the strange sea, as seen in Fig. 5b and 6a. Depending on details of the model, we get the 
range 0.0010 < S"" < 0.0023 for this asymmetry. 

This is especially interesting in connection to the NuTeV anomaly jTUI. NuTeV found, 
based on the observable = ^1!^I^n^I^-x)Z%''^nZI+^^^ = - f I = I - sm^^w, a 3a 
deviation of sin^ Ow compared to the Standard Model fit: sin^ ^^^."^"^ = 0.2277 ± 0.0016 
compared to sin^ = 0.2227±0.0004. However, an asymmetric strange sea would change 
their result, since u only have charged current interactions with s and u with s. Using the 
folding function provided by NuTeV to account for their analysis, the s-s asymmetry from 
our model gives a shift -0.0024 < Asin^^iy = dx s~{x)F{x) < -0.00097, i.e. the 
discrepancy with the Standard Model result is reduced to between 1.6cr and 2.4cr, leaving 
no strong hint of physics beyond the Standard Model. 

In our model, it is also natural to consider fluctuations involving heavy quarks. Assum- 
ing that the hadronic fluctuation description is still valid for the proton fluctuating into 
charmed baryon-meson pairs, the lightest such fluctuations are p AqD and p —* p J/ip. 
Taking these fluctuations into account implies that there should be an intrinsic charm 
component in the proton at intermediate x ~ 0.4. This component is quite different from 
the purely perturbative charm component from gluon splitting, which falls steeply with x 
(see Fig. 7). The intrinsic charm component should be present from the scale where the 
momentum transfer can realize the charmed fluctuation. 

The resulting distributions of c and c at the starting scale are shown in Fig. 6. In 
Fig. 6a the distributions from the \AcD) fluctuation are shown compared with the s and 
s distributions. The fluctuation parameters {aH and aq) are taken to be the same as for 
the light quark fluctuations, just as for the strange sea. However, the sensitivity of the 
result on the precise values of these parameters is small. The normalization is in Fig. 6 
taken to be oc 1/AMbm (to be discussed below), in order to easily compare the shapes 
of the strange and charmed sea. In this case, there is an asymmetry between the c and c 
distributions, similar to that of the strange sea, but much smaller due to the similarity in 
mass between the Ac (2285 MeV) and the D (1865 MeV). Fig. 6b shows the distributions 
from the \p J/ip) fluctuation. They are very similar to that of the \AcD) distributions, 
except that the asymmetry is missing since both the c and the c are here in the meson. 
From Fig. 6 it is clear that an investigation of the intrinsic charm distributions from our 
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FIGURE 7. Data for the charm structure function F2 from the EMC experiment [T2j, com- 
pared to our model results with the best-fit normalization (solid curve). The dotted curves show 
the perturbative QCD prediction for 7(7 — > cc from with three different choices for the renor- 
malization and factorization scales. For comparison, the intrinsic charm distribution of 11 is 
also shown, with the 0.7% normalization allowed according to |13j (dashed curve). 

model is not much affected by the precise nature of the dominating fluctuation mode, and 
in the following we will use only the \KcD) fluctuation. Note that the shape of our intrinsic 
charm component is somewhat different from that in the intrinsic charm model by Brodsky 
et al. jTTj, which is based on partonic fluctuations p — > uudcc (see Fig. 7). 

The normalization used in Fig. 6 corresponds to a relative importance between differ- 
ent mass states proportional to 1/AMbm, as suggested by the strange sea fit to CCFR 
data. This normalization would give an intrinsic charm component ((c + c)/2 integrated 
number density) of 0.9%. However, it might be more appropriate to use a normaliza- 
tion oc 1/ AM'^j^.j (as given by old-fashioned perturbation theory) compared to the strange 
fluctuations, corresponding to an intrinsic charm component of 0.18%. 

The only experimental data for the large-x charmed structure function F| = |x(c-|- c) 
(at leading order) comes from the EMC experiment 12 , which measured charmed hadron 
production in muon-proton scattering. There, an intriguing excess was found in the largest 
X bins, compared to the perturbative photon-gluon fusion expectation. A later analysis 
gave further evidence that the excess cannot easily be attributed to standard perturbative 
production channels JHl- Intrinsic charm was immediately suggested as an explanation 
for the excess, but the shape of the charm distribution in the original intrinsic charm 
model is not optimal to explain the EMC excess. In Fig. 7, the EMC data is shown 
in bins of u = Q'^/2MpX, together with the result from our model, evolved in using 
NLO QCD evolution jHj. Here we use the best-fit normalization, corresponding to 0.45% 
intrinsic charm. This lies between the two normalizations discussed above, which should 
not be surprising since the energy denominator only gives an order-of-magnitude estimate. 
For comparison, we also show the intrinsic charm distribution of ^T] with the largest 
normalization allowed by the EMC data according to 13j (0.7% intrinsic charm), and the 
perturbative photon-gluon fusion results from [TB] . 

As can be seen from Fig. 7, the shape of the intrinsic charm distribution in our model 
seems to fit the data very well, giving an enhancement at precisely the right values of x. 
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Unfortunately, the statistics of the EMC result is too small to allow any discrimination 
between different models for intrinsic charm, and measurements of the charm structure 
function at HERA (Hl^^ and ZEUSjini) values of x to contribute to 

our understanding of intrinsic charm. If a future experiment would measure the large-x 
charm structure function with large statistics, it would be very interesting to get a decisive 
verification of the presence of intrinsic charm in the proton. 

Acknowledgments: I would like to thank the organizers for the opportunity to talk 
at DIS'05, and Stan Brodsky and Gunnar Ingelman for very interesting discussions. 
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